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Abstract 

We show that fluctuations in elliptic anisotropy in heavy-ion collisions 
can be used to significantly reduce the contribution of transverse-momentum 
conservation, and of all background effects independent on the orientation 
of the reaction plane, from an observable of the chiral magnetic effect. We 
argue that for a given impact parameter, the magnetic field is approximately 
independent of the fluctuating shape of the fireball. 



1. The main feature of the chiral magnetic effect (CME) [T] is the existence of an 
electric current parallel (or antiparallel) to the direction of magnetic field produced 
in heavy- ion collisions [2], see also [I], [3|. Consequently, charge separation can be 
observed in the direction perpendicular to the reaction plane. Recently, evidence for 
the CME was found in lattice QCD calculations [4J. A detailed discussion of this 
effect was given in Ref. [5]. 

As discussed in Ref. [6], such charge separation can be observed via the following 
two-particle correlator 

7 = (cos(0i + 2 - 2*Rp)) , (1) 

where ^rp, 0i, and 2 , respectively, denote the azimuthal angles of the reaction 
plane, and the two charged particles produced. Alternative observables were pro- 
posed in Refs. (TJ [8] (see last section). Recently, the STAR collaboration measured 
7 and their result [HI HO] is qualitatively consistent with the CME expectation; how- 
ever, the interpretation of these data still are debatable, see e.g., [TT], [T51 dSl EEU 
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fT5l fT6l fT71 fl8] . The interpretation of experimental results is not simple because 
practically all two-particle correlations contribute to 7, as seen from (here ^pp = 0) 

7 = (cos(0i) cos(0 2 )) - (sin(0i) sin(0 2 )) • (2) 

Indeed, in the presence of elliptic anisotropy [Tj5] , both terms can differ even though 
the correlation mechanism itself is not directly sensitive to the orientation of the 
reaction plane. 

In this paper, we discuss this problem in detail. In the next Section, we undertake 
explicit calculations to quantify the contribution of elliptic anisotropy to 7. Next, we 
argue that in peripheral heavy-ion collisions (where the CME is considered to have a 
maximum strength [Tj) we expect large fluctuations in elliptic flow (t> 2 ) such that it 
allows us to deduct the v 2 -driven background from 7. We also show that the changing 
shape of the fireball, at a given impact parameter, does not change the contribution 
of the CME to 7. Finally, we present a Monte Carlo model, wherein we evaluate 
the contribution of transverse-momentum conservation to 7 at vanishing elliptic 
anisotropy (f 2 — > 0). In the last section we give our comments and conclusions. 

2. By definition 

_ / P2(01, 02, Si, x 2 , ^Ipp) cos(0i + 2 - 2# RP )d(j) 1 d(j)2dx 1 dx2 ^ 

/ P2(</>1, 02, X 2 , ^ R P )d(j) 1 d(j)2dx 1 dX2 

where p 2 is the two-particle distribution at a given angle of the reaction plane, 
typp. To make our notation shorter, we denote x = (pt, rj) and dx = ptdptdrj, where 
p t is the absolute value of transverse-momentum, while 77 is pseudorapidity. The 
distribution p 2 can be expressed via the correlation function C 

p 2 (01,0 2 ,Xi,X 2 ,^ fiP ) = p((p 1 ,X 1 ,^ R p)p((f)2,X2,^Rp)[l + C(0i,0 2 ,Xi,X 2 )], (4) 

with the single-particle distributiorJl] 

p{<f>, x, Vrp) = + 2v 2 {x) cos (20 - 2ttju0], (5) 

where po{x) does not depend on and ^ pp. We study only those correlations that 
do not depend on the reaction plane, i.e., C only depends on A0 = 0i — 2 . Next 
we expand C in a Fourier series 

Eoo 
a n (xi,a; 2 )cos(nA0) , (6) 
n=0 

where a n (xi,x 2 ) does not depend on 0i and 2 . Substituting ([6]) and (TJJ into Eq. 
Q, we obtain 



7 = 2iV2 ' Po( x ^Po( x 2)ai(x 1 ,X2)[v 2 (x 1 ) + V2(x 2 )]dx 1 dx 2 , (7) 



1 Higher v n results in terms proportional to f2f4, V4VQ, etc., and can be neglected, see Eq. ([7|). 
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where N — J po(x)dx and we assume that 1 + a n ~ 1. As seen from Eq. (JTj), all 
correlations with non-zero ai(xi,x 2 ) contribute to 7, even if the underlying corre- 
lation mechanisms do not depend on the orientation of the reaction plane. This 
finding explains why transverse- momentum conservation [Tlj [HJ IE] , local charge- 
conservation [12], resonance- (cluster-) decay [15], and all other correlations with A0 
dependence contribute to 7. However, as pointed out recently in Ref. [2D], those 
correlations can be removed from 7 by taking only those events where V2(x) ~ 0. We 
note that v 2 (x) is defined solely through Eq. fl5]), and it can be positive or negative. 
Taking the CME into account^ 

p x (fa x, V RP ) = ^ [1 + 2v 2 (x) cos (20 - 2V RP ) + 2 X d(x) sin (0 - * RP )] , (8) 

Z7T 



we obtain, 
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~1 2 



d(x)p (x)dx 



+ 



2iV 2 



Po(^i)po(^2)ai(a;i, x 2 )[v 2 (xi) + v 2 (x 2 )]dx 1 dx 2 , (9) 



wherein the first term represents the CME, and the second term is the elliptic- 
anisotropy-driven background. 

3. We expected that elliptic anisotropy would be correlated with the participant 
eccentricity e 2 [2T]. In the center of mass of the wounded nucleons, e 2 is given by 



'(E i r l 2 cos(20 l ))^ + (E l r J 2 sin(20 l ))^ 
e 2 = ^-2 > ( 10 ) 

wherein the wounded nucleons are characterized by their radii, r,, and their az- 
imuthal angles, fa. 

In Fig. [H we present the calculate^ e 2 distribution in Au+Au collisions at 
y/s = 200 GeV with b = 10 fm (the impact parameter). Accordingly, even at b = 10 
fm, we obtain a broad range of e 2 (and v 2 j& that can be used to significantly change 
the background present in the second term of Eq. (j9j) . 

However, we wanted to remove this background under the condition that the 
contribution from the CME is approximately unchanged. To verify this, we calcu- 
latecfl the out-of-plane component of the magnetic field By at t — (time) as a 



2 The value of \ flips between —1 and +1 so that | J2 X Px = P> defined in Eq. (JSJ). 

3 We use the Monte Carlo Glauber calculation with standard parameters [22] . 

4 In 3D event-by-event hydrodynamics [23], (v 2 ) = 0.08 and - (w 2 ) 2 ] 1/2 = 0.04 for 40-50% 



centrality (b w 10 fm) in Au+Au collisions. 

5 The details of our calculation are presented in Ref. 
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Figure 1: Normalized density distribution of the participant eccentricity 62 at the 
impact parameter b = 10 fm in Au+Au collisions at y^i = 200 GeV. 



function of €2- As seen in Fig. [2j the magnetic fields from wounded- and spectator- 
protoni are approximately constant (in comparison to v 2 that scales linearly with 
62) in the broad range of €2- We conclude that fluctuating in peripheral collisions 
will allow us to study the V2 dependence of 7 at an approximately constant strength 
of the CME. 
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Figure 2: The out-of-plane component of magnetic field at b = 10 fm in Au+Au 
collisions at a/s = 200 GeV as a function of the participant eccentricity €2- Contri- 
butions from the wounded- and spectator- protons are depicted separately. 

In Figs. [1] and [2] we made our calculations at a given impact parameter. In an 
actual experiment, we can select our peripheral collisions, e.g., by the number of 

6 Both sources are expected to contribute differently during the evolution of the fireball [T|. 
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particles produced at midrapidity. We consider that our conclusions also will hold 
in this situation. 

4. As an example, we calculated 7 as a function of v 2 in a model with only 
transverse-momentum conservation and elliptic anisotropyjj As shown in [TTj [HI 
OS], this effect contributes significantly to 7 and reasonably describes the p t and 77 
dependence. 

We sampled N a u = 50 particles with p t according to the thermal distribution 
p t e~ Pt/T with 2T = 0.45 MeV, and according to 1 + 2v 2 (p t ) cos(20). We took 
v2i.Pt) = 0.14p t for p t < 2 GeV and v 2 {p t ) = 0.28 for p t > 2 GeV, so that we obtained 
the integrated v 2 ~ 0.06. Next, we imposed transverse-momentum conservation^ 
and calculated 7 as a function of selected v 2 = j^— cos(20j), where N obs is the 

number of observed particles (selected randomly from N a u). Owing to the statistical 
fluctuations, we obtained a broad range of v 2 that allows us to test Eq. (j7|). 




transverse momentum conservation 



0.00 ^05 OlO^ 0.15 0.20 

v 2 



Figure 3: The contribution of transverse- momentum conservation to 7 as a function 
of v 2 . The momentum is conserved for N a u particles and 7 is calculated for N b s 
(observed) particles, selected randomly from A^. 

As seen from Fig. |3l 7 — > as v 2 — > for N a u = N Q b s = 50. However, it might be 
surprising that for N b s = 20, 7 becomes slightly positive as v 2 — > 0. In Ref. p4] the 
contribution of conserving transverse-momentum was calculated analytically, and 
was shown to scale not only with v 2 in the region where particles are observed, but 

7 The calculations presented in this section are for illustrative purposes only, and should not be 
compared with the STAR data. 

8 We accepted only those events where the total (summed over N a u particles) \pt,x\ < 0-5 GeV 
and \p t ,y\ < 0.5 GeV. 
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also with t>2 in the whole phase-space 



1 



[2v 2 ,n - v 2 ,f] , 



(11) 
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where v 2 ~ / po(x)v2(x)ptdx, v 2 ~ / po(x)v 2 (x)pjdx and x = (p t , 77), c?x = p t dp t dr]. 
The indexes F and f2 indicate that integrations are performed over full phase-space 
(F) or the phase space wherein particles are measured (O), respectively. When we 
calculate 7 for all produced particles, Q = F, then v 2 ^ = implies v 2j n(p t ,ri) = 
and t^n = v 2) p = 0. Consequently, 7 = at u 2 ,n — 0. However, if we measure only 
a fraction of all particles, then for v 2) q = (and v 2} q = 0), v 2 ,f can differ from zero 
(positive) and 7 ~ v 2 ^/N a u > 0, as seen from Eq. ( ITT]) . 

5. Several comments are warranted: 

(i) Very recently, an observable related to 7 was studied as a function of elliptic 
anisotropy [20] . This finding argued that for mid-peripheral Au+Au collisions 7 < 
at v 2 = for same-charge pairs, which is consistent with the CME. 

(ii) Even if 7 < at v 2 = for same-charge pairs, it does not imply the existence 
of the CME. It only indicates the presence of some correlation mechanism that 
explicitly depends on the orientation of the reaction plane. To measure the CME, a 
different observable is needed, e.g., the multiparticle charge-sensitive correlator [7J 
or direct measurements of the electric dipole [8]. 



(iii) As argued in Ref. [22], central U+U collisions also can be used to distinguish 
between effects driven by elliptic anisotropy and the CME. In the present paper, we 
considered only peripheral collisions, where the CME and fluctuations in v 2 are 
expected to be the most visible. Both methods can be used independently to reduce 
the contribution of elliptic anisotropy. 

(iv) In this paper, we proposed a way to remove the elliptic-flow-induced back- 
ground from the correlator 7. However, we note that in Ref. [7] a new multiparticle 
charge-sensitive correlator C c was proposed that is insensitive to correlations due to 
(elliptic) flow, jets, or momentum conservation. The measurement of C c together 
with 7 (with removed background) could provide important information about a 
possible signal of local parity violation. 

In summary, we demonstrated that fluctuations in elliptic anisotropy in pe- 
ripheral heavy-ion collisions can be used effectively to reduce the contribution of 
U2-induced correlations from the two-particle correlator ([T|). We showed that at a 
given impact parameter, the magnetic field produced in heavy-ion collisions depends 
weakly on the participant eccentricity, in contrast to the value of v 2 . We also dis- 
cussed the contribution of transverse-momentum conservation to 7 at a vanishing 
v 2 . Preliminary experimental analysis [20] suggests the presence of a correlation 
mechanism that does not scale with elliptic anisotropy. 
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